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(54) Interpolation circuit fdr an Incremental encoder 



(57) Output of an encoder are sampled by A/D con- 
verters 11a and 1 lb to be converted to N bits A-phase 
and B-ptiase digital data DA and DB. In a look-up table 
memory 12, reference pfiase angle data of phase divi- 
sions and average gradient vectors of cfianges in phase 
angle data within the phase divisions are stored, the 
phase divisions being addressed by the high order NU 



bits of data DA and DB. An arithmetic circuit 13 deter- 
mines a vector inner product of an average gradient 
vector and phase-interpolating data represented by the 
low order NL bits of the data DA and DB to add the 
resultant to a phase angle data, thereby outputting an 
interpolated phase angle data. 
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Description 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

This invention relates to an interpolation circuit performing digital interpolation processing to two-phase sinusoidal 
signals of an encoder, which detects position, angle, velocity, angular velocity and so on, to ot>tain phase angle data 

with high resolution. ' 

10 

Related Art 

Since there is a machining limit in grid spaces formed in a scale of the encoder, a space period of the sinusoidal 
signals from the encoder should be subdivided to be interpolated in order to measure a space finer than the scale grid. 
IS Therefore, various interpolation circuits have been used previously. 

Fig.4 is an example of the interpolation circuit by previous digital processing. Two-phase sinusoidal signals INA and 
INB of a phase A and a phase B shifted by 90° each other are output from an encoder 1 . and then sampled with a pre- 
determined fl-equency by A/D converters 2a and 2b to be converted to digital data DA and DB respectively. Phase angle 
data at each sampling point are prepared and stored in a look-up table memory 3 in advance using an arc tangent func- 
20 tion(ATAN) based on the following formula. 

u = ATAN (DB/DA) (1) 

Therefore, the phase angle data u at each sampling point is obtained by reading the look-up table memory 3, defin- 
25 ing ttie digital data DA and DB as addresses x and y respectively Furthermore, the phase angle data in is entered in a 
two-phase square wave data generation circuit 4, whereby digital two-phase square wave data OUTA and OUTB are 
obtained. 

When Intending to obtain a sufficient interpolation number in the interpolation circuit in which the look-up table 
memory is used as desaibed above, capacity of the look-up table memory becomes extremely large. For example, 
30 when the digital data DA and DB obtained by A/D converters 2a and 2b are N bits, the size of the address space spec- 
ified by the digital data DA and DB is N x N. When a required interpolation number is defined as I, the Integer J beyond 
log2l is required to the phase angle data. In such a condition, the memory capacity required in the look-up table memory 
becomes 22"^ X J. 

For example, the memory configuration where tiie phase angle data obtained by dividing one period into 400 are 
35 addressed with the address of 8 bits x 8 bits becomes as Fig.5. When a data length is defined as J = 9, the memory 
capacity becomes 2^ x 2^ x 9 = 589,824 bits. 

As a method for reducing the capacity of the look-up table memory in the encoder interpolation circuit of a digital 
system as described above, a method for reducing the address space have been proposed, paying attention to the sym- 
metry that tiie phase angle data is repeated periodically wifliin an address space with exception of a sign (Japanese 
40 Laid-open Patent Publication NaHei.3-68812). 

SUMMARY OF THE INVENTION 

The object of this invention is to provide an interpolation circuit of an encoder wherein the phase angle data defined 
45 by the high order bits of the digital data obtained by A/D conversion is designed to be interpolated utilizing the low order 
bits thereof whereby the capacity of ttie look-up table memory is able to be reduced. 

The interpolation circuit of the encoder according to the present invention comprises: an A/D conversion circuit for 
sampling A-phase and B-phase sinusiodal signals, which are output from an encoder as phase-shifted by 90° each 
ottier, witii a given frequency to convert ttie resultant to A-phase and B-phase digital data witti N bits, respectively: a 
so look-up table memory which stores reference phase angle data con-esponding to a plurality of phase divisions and inter- 
polation factors prepared for every phase divisions for interpolating the reference phase data to obtain an interpolated 
phase angle data, the phase divisions being obtained by dividing an phase space, two-dimensional address of which is 
assigned by the A-phase and B-phase digital data, to have two<limensional addresses specified by ttie high order NU 
bits data in ttie A-phase and B-phase digital data, and in which ttie high order NU bits data is entered as ttie two-dimen- 
55 sional address signals; and an ariflimetic circuit for obtaining ttie interpolated phase angle data by processing ttie ref- 
erence phase angle data and the respective interpolation factor read out from ttie look-up table memory and a phase- 
interpolating data represented by the low order NL (=N-NU) bits data in ttie A-phase and B-phase digital data. 

In the present invention, the phase space usually required to a conventional look-up table memory corresponding 
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to a required interpolation number is divided into a plurality of phase divisions. In the look-up table memory in the 
present invention, the phase angle data con-esponding to the number of the phase divisions respectively defined by the 
high order NU bits of the A-phase and B-phase digital data are stored, and the reference phase angle data Is read to 
he interpolated by the low order NL bits sata. The Interpolation factors used in the interpolation of the phase angle data 

5 are stored in the look-up table memory as well as the phase angle data defined by the high order NU bits data. There- 
fore, the interpolated phase angle data can be obtained by simple interpolation operation in such a manner that the 
required interpolaton number is maintained, while the capacity of the look-up table memory is reduced by a large 
amount as compared to the case that the look-up table memory' has an full address of N x N and the entire phase angle 
data defined by this address are stored. 

10 Concretely, since the two-dimensional address is used Ibr processing the A-phase and B-phase sinusiodal signals, 
and the phase-interpolating data and the interpolation factor are expressed as two-dimension data, vector operation is 
used for interpolation operation. For example, in the present invention, an average gradient vector of a change in phase 
angle data within each small address space as the interpolation factors as well as the reference phase angle data of 
the reference position of each phase divisions being addressed by the high order NU bits of the A-phase and B-phase 

15 digital data are stored in the look-up table memory. In this case, such an interpolation operation is performed that a vec- 
tor inner product of the phase-interpolating data composed of the low order NL bits of the A-phase and B-phase digital 
data and the average gradient vector read from the look-up table memory is produced, and then the resultant is added 
to the reference phase angle data read from the look-up table memory. As a result, a deviation (an error) of the phase 
angle data Interpolated according to the present invention from the phase angle data which is obtained when using the 

20 look-up table memory having the full address becomes sufficiently small. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The embodiments according to the present invention will be described refening to the following drawings. 
25 Fig. 1 shows an interpolation circuit of an encoder according to an embodiment of the present invention. 

Fig.2 is a view for illustrating principles of interpolation of phase angle data according to tiie embodiment. 

Fig.3 is a view for illustrating a relationship between a reference phase angle data and an average gradient vector 
according to the embodiment. 

Fig.4 shows an interpolation circuit of a prior encoder. 
30 Fig.5 shows an address and a phase space of a look-up table memory in the interpolation circuit in Fig.4. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Fig.1 shows an interpolation circuit of an encoder according to an embodiment of tiie present invention. Although 

35 it is no matter what principle is used for the encoder to which the interpolation circuit is applied, an optical encoder and 
a magnetic encoder using MR elements are preferably used as the encoder. The A-phase and B-phase sinusoidal sig- 
nals INA and INB output from the encoder are sampled by A/D converters 1 1 a and 1 lb witti a predetermined frequency, 
respectively, whereby tiie A-phase and B-phase digital data DA and DB are obtained. A look-up table memory 12 con- 
stituted by nonvolatile memory such as ROMs is used In order to determine, phase angle. The look-up table memory 

40 12 is accesses by the digital data DA and DB. This basic configuration is similar to prior art. However, in tills embodi- 
ment, the phase angle data, number of which is determined by use of predetermined high order bits of the data DA and 
DB as two-dimensional address signals, are stored in the look-up table memory 12. 

Namely, in the case of the embodiment, the data DA and DB are N bits (N = 8 in Fig.1 ), respectively The high order 
NU bits (NU = 5 in Fig.1) of these data DA and DB are entered into tiie look-up table memory 12 as two-dimensional 

45 address signals. In the memory 1 2, tiie reference phase angle data uO con-esponding to the respective reference posi- 
tions (xO, yO) of the respective phase divisions "a" are stored witii interpolation factors which are prepared con'espond- 
ing to the respective reference phase angle data for phase angle interpolating. The phase divisions "a" are obtained by 
dividing an phase space, each two-dimensional address of which is assigned by tiie A-phase and B-phase digital data 
to have two-dimensional addresses specified by ttie high order NU bits data in the A-phase and B-phase digital data. 

50 The size of each phase division "a" is NL x NL required to store the phase angle data of 2^'^'- words (64 words in the 
case of NL = 3) in tiie case of using the digital data DA and DB of N bits as the full address, as shown in Fig.3. 

As is described below, average gradient vector data of tiie changes in ttie phase angle data wittiin ttie respective 
phase divisions "a" are used as ttie interpolation factors being stored as well as ttie reference phase angle data uO. It 
is appreciated that each of the reference positions of ttie reference phase angle data uO is a predetermined position in 

55 the respective phase division "a". 

Data Ax and Ay of tiie low order NL bits (NL = 3 in Fig.1) of ttie respective data DA and DB are entered into an arith- 
metic circuit 13 for performing interpolation operation as phase-interpolating data, respectively In the aritiimetic circuit 
13, operation for interpolating ttie phase angle data is performed based on ttie reference phase angle data and ttie 
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interpolation factors (i.e., the average gradient vectors) read out from the look-up table memory 12 and the phase-inter- 
polating data con-esponding to the low order bits data Ax and Ay. 

Referring to Fig.2, the principle of interpolation of phase angle data according to this embodiment will be described. 
Fig.2 shows only the upper right quadrant of the phase space wherein the data DA and DB are represented as 

5 addresses x and y. In Fig.2, uO is a certain reference phase angle data read out from the look-up table memory 12. and 
A u is a vector of the phase-interpolating data defined by the data of NL bits x NL bits that is to be included in the 
phase division "a ". A phase angle data u to be obtained in the case of using the full address of N bits x N bits is obtained 
by interpolating the reference phase angle data uO using the vector A u of the phase-interpolating data. The amount 
of changes (the gradient) in the phase angle data within each phase division "a" is determined by the reference position 

10 (xO, yO) of the phase divisions "a", and is substantially uniform. Therefore, as described above, the average gradient 
vector Ki of the phase angle data in the phase divisions "a" is determined in advance to be stored as the interpolation 
factor In the look-up table memory 12 as well as the reference phase angle data uO. An interpolated phase angle data 
u can be obtained by performing a vector inner product operation of the vector A u of the phase-interpolating data and 
the average gradient vector k: , and by adding the result to the reference phase angle data uO, as represented by ttie 

75 following formula (2). 



u = uO+l<: xAu (2) 

20 

X-, y-components of the phase-interpolating data vector A u are Ax, Ay. respectively. When x-, y-componenis of 
the average gradient vector icare defined as kx. ky. respectively, the formula (2) becomes the following formula (3). 

25 u = uO + (kx- Ax + ky Ay) (3) 

The interpolation operation using the vector inner product described above is performed in the arithmetic circuit 13 
In Fig. 1 . A first and a second multipliers 1 3 1 and 1 32 multiply the components kx and ky of the average gradient vector 
K , which are read out from the memory 12, and components Ax and Ay (i.e., components of the phase-interpolating 
30 data) represented by the low order NL bits of the digital data DA and DB, respectively. A first adder 133 adds the multi- 
plied results of the multipliers 131 and 132. The phase-interpolating data vector A u can be obtained by the adder 1 33. 
A second adder 134 adds the phase-interpolating data vector A u to the reference phase angle data uO read out from 
the memory 12. As a result an Interpolated phase angle data u that is to be required can be obtained. 

The interpolated phase data u Is. tor example, processed In a two-phase square wave data generation circuit 1 4 In 
35 a conventional manner, whereby two-phase square wave data OUTA and OUTB are output. 

For example, the components kx and ky of the average gradient vector K in the phase division "a" described above 
are given by a formula (4) described below, defining a center of the phase division "a" as a typical point (x1, y1) as 
shown in Figure 3. 

40 kx = -A.y1/(x1%y1^) (4) 

ky = +A'y1/(x1^+y1*) 

In formula (4), when the number of Interpolation is defined as I, "A" is; A = 1/2 n .The formula (4) is introduced 
45 using a mathematical formula of a scalar field. That is. when data represented by the scalar field having a certain spatial 
distribution of the changes is defined as u = A • (|i , and x-, y-components of a change portion vector A u of the data 
arc defined as Ax. Ay. respectively, "u" can be given by the following partial differential equation (5). 

u = uO(xO,yO) + (5u/5x)Ax + (5u/5)0Ay (S) 
^ = uO (xO. yO) + (- Ay1)/(x1 % y1 ^) + (+ Ax1)/(x1 ^ + y1 ^) 

Based on the mathematical formula, then the data to be stored in the memory is defined as A = 1^ it (I is defined 
as the number of interpolation), the components of the average gradient vector of the formula (4) are obtained. 
55 As is described in Fig. 3, it is preferable that the typical point (x1 , y1) determining the average gradient vector K is 
defined as a center of the phase division "a" In order to improve the accuracy of approximation. However, since the 
changes within the phase division are small, a large approximation error can not be occurred even when other typical 
poirrts are selected. For example, the reference position (xO, yO) can be selected as a typical point. 
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Moreover, as the other simple method determining the average gradient vector of the changes in the phase angle 
data, for e)rample, the average gradient vector can be obtained by calculating a difference between two typical points of 
adjacent phase divisions as follows. 

la = {u(xO + 2NL,yO)-u(x0.yO)}/2^'- (6) 
ky = {u (xO. yO + 2NL) - u (xO, y0)}/2 

It will be described that the approximation error of interpolation of the phase angle data using the aforementioned 
average gradient vector k: is small, concretely providing the examples of numeral calculations. 

When the number of interpolation I is 400, the phase angle data u at several positions y/x (= DB/DA) in a range in 
which the signal amplitude becomes approximately one-half of the minimum, is determined as described In the follow- 
ing Table 1. 



. TABLE 1 

Phase angle data u 



19 


18.37204 


18.1047 


17.84443 


17.59154 


18 


17.45404 


17.19849 


16.95013 


16.70868 


17 


16.52854 


16.28528 


26.04893 


15.8192 


16 


15.59583 


15,36517 


1 5.1411 


14.92337 


y/x 


64 


65 


66 


67 



In TABLE 1 , when the difference data A u is determined, defining (x, y) = (64, 1 6) as the reference position, resulting 
in as Table 2. 



TABLE 2 

Difference data A u 



19 


2.776216 


2.508671 


2.248604 


1.995715 


18 


1.858216 


• 1.602661 


1.354306 


1.112858 


17 


0.932709 


0.689453 


0.453103 


0.223374 


16 


0 


-0.23066 


-0.45472 


-0.67246 


y/x 


64 


65 


66 


67 



Next, when partial differential quotients 5 u/ S X and S u/ 5 y of each phase angle (these quotients correspond to ttie 
gradient vector components kx and ky) are determined firom tiie phase angle data u of the TABLE 1, resulting in as 
TABLE 3 and Table 4 described below. 
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d w s X 





19 


-0.27139 


-0.26375 


-0.25643 


-0.24940 




18 


-0.25926 


-0.25190 


-0.24485 


-0.23809 


10 


17 


-0.24681 


-0.23975 


-0.23299 


-0.22651 




16 


-0.23405 . 


-0.22731 


-0.22086 


-0.21467 


15 


y/x 


64 


65 


66 


67 



TABLE 4 



8 u/ S y 





19 


0.91415 


0.90232 


0.89075 


0.87945 




18 


0.92180 


0.90966 


0.89780 


0.88621 


SB 


17 


0.92916 


0.91671 


0.90456 


0.89271 




16 


0.93621 


0.92346 


0.91103 


0.89892 




y/x 


64 


65 


66 


67 



The interpolated phase angle data u = uO + kx* Ax + lv Ay is determined by using the reference phase angle 
35 uO of y/x = 16/64, using the partial differential quotients 5 u/ 5 x = - 0.24238 and 5 u/ 5 y = 0.907181 in x1 = 65.5, y1 = 
17.5 as the average gradient vector components kx and 1^ from TABLE 3 and TABLE 4, resulting in as TABLE 5 
described below. 

TABLES 

40 



u = uO + kx- A X + ky- A y 





19 


18.31737 


18.07499 


17.83262 


17.59024 


45 


18 


17.41019 


17.16781 


16.92544 


16.68306 




17 


16.50301 


16.26063 


16.01825 


15.77588 


SO 


16 


15.59583 


15.35345 


15.11107 


14.8687 




y/x 


64 


65 


66 


67 



55 

The results of ttie TABLE 5 correspond to the phase angle data being approximately obtained by ttie interpolation 
circuit of Vns embodiment. The errors between TABLE 1 and TABLE 5 are summarized as described In tiie following 
TABLE 6. 
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TABLE 6 



Error 



19 


-0.05467 


-0.02950 


-0.01181 


-0,00130 


18 


-0.04385 


-0.03067 


-0.02470 


-0.02563 


17 


-0.02553 


-0.02467 


-0.03067 


-0.04332 


16 


-0.00000 


-0.01172 


-0.03003 


-0.05467 


y/x 


64 


65 


66 


67 



From the numeral data described above, the errors of the approximate value obtained by the interpolation circuit of 
20 the embodiment are such small as less than ± 1 LSB of the number of interpolation 400. 

As is described above, when the phase angle data of the number of interpolation 400 was assigned in the lool<-up 
table memory being addressed by the address of 8 bits x 8 bits, the memory capacity was 2^ x 2^ x 9 = 589,824, defin- 
ing a data length of the phase angle data as 9. In contrast with this, in case of the embodiment, assuming that the ref- 
erence phase angle with the same data length and the gradient vector components kx and ky of 5 bits respectively are 
25 Stored in the look-up table designated by the high order bit of 5 bits x 5 bits, the memory capacity becomes 2^ x 2^ x (9 
+ 5 + 5) = 19,456 bits. Namely, the memory capacity can be compressed to one thirtieth in comparison with the above- 
described example. 

Such Nonvolatile memories as ROM, EPROM and EEPROM, and PLD, RAM or the like can be used as the look- 
up table memory 12 shown in Flg.1 . Moreover, the arithmetic circuit 13 performing interpolation operation can be real- 
30 ized tiirough software by a microprocessor or by a DSP except the digital circuit. 

According to the present invention, in tiie look-up table memory for interpolating, the phase angle data designated 
by the high order NU bits of the digital data obtained by converting output of the encoder by A/D converter and the aver- 
age gradient vector of the changes In the phase angle data determined in advance in response to the phase position of 
the reference phase angle data are stored. The interpolation operation is performed based on the data of the low order 
SB NL bits, and the reference phase angle data and the average gradient vector read from tiie look-up table memory, 
whereby the capacity of the look-up table memory can be largely compressed. 

Although the present invention has been shown and described with respect to best mode embodiments tiiereof, it 
should be understood by those skilled in the art that the foregoing and various other changes, omissions, and additions 
in the form and detail thereof may be made therein witiiout departing from ttie spirit and scope of tiie present invention. 
40 The entire disclosure of Japanese Patent Application No. 9-1111 23 filed on April 28, 1 997 including specification, 
claims, drawings and abstract are incorporated herein by reference in its entirety. 

Claims 

45 1. An interpolation circuit for an encoder, comprising: 

an A/D conversion circuit for sampling A-phase and B-phase sinusiodal signals, which are output from an 
encoder as phase-shifted by 90° each ottier, with a given frequency to convert ttie resultant to A-phase and B- 
phase digital data witt) N bits, respectively: 

50 a look-up table memory which stores reference phase angle data con-esponding to a plurality of phase divi- 

sions and interpolation Actors prepared for every phase divisions for interpolating the reference phase data to 
obtain an interpolated phase angle data, the phase divisions being obtained by dividing an phase space, two- 
dimensional address of which is assigned by ttie A-phase and B-phase digital data, to have two^limensional 
addresses specified by ttie high order NU bils data in ttie A-phase and B-phase digital data, and in which ttie 

55 high order NU bits data Is entered as the two-dimensional address signals: and 

an arithmetic circuit for obtaining the interpolated phase angle data by processing the reference phase angle 
data and the respective interpolation factor read out from the look-up table memory and a phase-interpolating 
data represented by the low order NL (=N-NU) bits data in the A-phase and B-phase digital data. 
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2. The interpolation circuit according to claim 1 , wherein the interpolation Actors stored in the look-up table memory 
are average gradient vectors of changes of the phase angle data predetermined in response to phase positions of 
the phase divisions. 

5 3. The interpolation circuit according to claim 2, wherein the arithmetic circuit comprises: 

mult'plicat'on means for determining a vector inner product of the phase-interpolating data and the average 
gradient vector read out fi-om the look-up table memory; and 

addition means for adding the vector inner product determined by the multiplication means to the reference 
10 phase angle data read out from the look-up table memory to output the interpolated phase angle data. 

4. The interpolation circuit according to claim 2, wherein the arithmetic circuit comprises: 

first and second multiplication means for performing multiplication of kx • Ax and ky • Ay, respectively, where 
IS Ax , Ay are x-, y-opponents of ttie phase-interpolating data, respectively, and kx, ky are x-, y-components of the 

average gradient vector read out from the look-up table memory, respectively; 
first addition means for adding ttie results obtained by ttie first and second multiplication means; and 
second addition means for adding the result added by the first addition means to the reference phase angle 
data read out from the look-up table memory. 

20 

5. The Interpolation circuit according to claim 2, wherein x-, y-components of the average gradient vector are obtained 
as partial differential quotients of the phase angle data at a center point within the respective phase division. 

6. The interpolation circuit according to claim 2, wherein x-, y-components of the average gradient vector are obtained 
25 as differentiated values between two typical points in ttie respective phase division and anottier phase division 

adjacent thereto. 
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FIG. 2 
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